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Abstract: In this work, we have continued the study of the Hawking radiation on 
the brane from a higher-dimensional rotating black hole by investigating the emission of 
fermionic modes. A comprehensive analysis is performed that leads to the particle, power 
and angular momentum emission rates, and sheds light on their dependence on fundamen- 
tal parameters of the theory, such as the spacetime dimension and angular momentum of 
the black hole. In addition, the angular distribution of the emitted modes, in terms of the 
number of particles and energy, is thoroughly studied. Our results are valid for arbitrary 
values of the energy of the emitted particles, dimension of spacetime and angular momen- 
tum of the black hole, and complement previous results on the emission of brane-localised 
scalars and gauge bosons. 
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1. Introduction 



The formulation, a few years ago, of theories JT], || postulating the existence of additional 
spacelike dimensions in nature have radically changed the landscape in gravitational physics 
(for some early works towards this direction, see Q). Driven by the motivation to address 
the hierarchy problem, these theories leave unchanged the usual 4-dimensional particle 
physics by localising the Standard Model (SM) particles on a 4-dimensional brane while 
allowing the gravitons to propagate in the bulk - the whole of the (4 + n)-dimensional 
spacetime. In the case of the scenario of Large Extra Dimensions jjj, in which the addi- 
tional spacelike dimensions are compact with size L, the latter feature introduces a new 
fundamental higher-dimensional scale for gravity M*, that can be significantly lower than 
the 4-dimensional scale Mp, thus increasing the fundamental Newton's constant - and the 
strength of gravitational interactions - by orders of magnitude. 

An immediate consequence of the above assumption is that particles with energies 
larger than M* can probe physics beyond the scale of quantum gravity while their collisions 
may trigger the appearance of strong gravitational phenomena and lead to the creation of 
extended objects (p-branes, string balls, string states, etc), rather than ordinary particles, 
including tiny black holes ||. This type of trans-planckian particle collision could be 
achieved in various environments, such as at ground-based colliders H or in high energy 
cosmic-ray interactions in the atmosphere of the Earth || (for an extensive discussion of 
the phenomenological implications and additional references, see the reviews [Q, g[ [|). The 
black holes produced with horizon radius <C L are submerged into the higher-dimensional 
spacetime, which modifies their properties (horizon radius, temperature, entropy, life-time 
etc.) [[R]]. Nevertheless, they are still expected to go through the same stages in their life 
as their 4-dimensional counterparts: a short balding phase, during which the black hole will 
shed all additional quantum numbers apart from its mass, angular momentum and charge, 
the more familiar spin-down phase, during which the black hole will lose mainly its angular 
momentum, then the Schwarzschild phase, when the black hole gradually loses its actual 
mass, and finally the Planck phase, where the black hole has reduced to a quantum object 
with unknown properties. 

During the two intermediate phases, the spin-down and Schwarzschild, the black hole 
loses energy through the emission of Hawking radiation [jllj (and through superradiance 
during the former phase). This will take place via the emission of elementary particles 
both in the bulk and on the brane, and it will be characterised by a very distinct thermal 
spectrum. It is this emission spectrum that upon successful detection will provide the 
most convincing evidence not only for the creation of microscopic black holes but of the 
very existence of extra spacelike dimensions. For this reason, the emission spectra from 
a higher-dimensional black hole have been under intensive study over the last few years. 
Having the simplest gravitational line-element, the Schwarzschild phase was the first one 
to undergo a systematic study both analytically 12, 13, [l4[] and numerically [15]. Both 
bulk (scalars) and brane (scalars, fermions and gauge bosons) channels were studied, and 
the derived emission spectra were shown to depend strongly on the number of additional 
spacelike dimensions in nature, both in terms of the amount of energy and number of 
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particles produced per unit time and of the type of particles produced. More general 
spherically-symmetric backgrounds were also studied, such as Schwarzschild-de Sitter [16], 
Schwarzschild-Gauss-Bonnet jT^] and Reissner-Nordstrom [18] backgrounds, as well as the 



effect of the mass of the emitted particles on the spectrum [19]. Recently, the graviton 
emission rates by a higher-dimensional Schwarzschild black hole were also investigated by 
a number of authors p(| . 

After the spherically-symmetric Schwarzschild phase was almost exhaustively studied, 
the investigation of the form of the radiation spectra for a higher-dimensional rotating 
black hole was initiated. Two early analytical studies 22] were supplemented, during 
the last year or so, by a number of works attempting to shed light on this more technically 
involved case. Exact numerical results for the emission of Hawking radiation from a (4+re)- 
dimensional rotating black hole, in the form of scalar fields on the brane, were first presented 
in [23]. In that work, the exact value of the angular eigenvalue, that couples the angular and 
radial part of the equation of motion of the scalar field and does not exist in closed form, 
was found numerically, and the radial equation was also solved by employing numerical 
analysis. A comprehensive study of the fluxes of particles, energy and angular momentum 
for the emission of scalar fields on the brane soon followed ]^4|] , that once again presented 
exact results for the emission rates including the angular distribution of particles. A similar 
extensive analysis, focussed on the emission of gauge bosons on the brane by a rotating 
black hole, appeared in [25|. In both of the latter works, it was shown that all emission rates 
for scalars and gauge fields are greatly enhanced by both the dimension of spacetime and 
the angular momentum of the black hole. Results produced independently in |2^, |27j, for 
these two species of particles, reveal a good agreement between the two sets of analyses. 
Additional studies in the literature featured the study of the superradiance effect for a 



higher-dimensional black hole [28] and of the emission in the bulk |2£ 

Recently, and while this work was in progress, a study |27j] produced power and angular 
momentum emission spectra from a rotating black hole for the last species of brane-localised 
particles, i.e. fermions. The radial part of the spinor equation was solved numerically 
with the angular eigenvalue given by an analytical expansion (with terms kept up to 6th 
order). The power and angular momentum emission spectra were found to exhibit the 
same characteristics as the ones for scalars and gauge bosons. The present work aims to 
complete the above analysis by presenting all emission spectra (particle, energy and angular 
momentum) as well as total emissivities (emission rates integrated over the frequency) and 
angular distribution of the emitted radiation. To this end, both the radial and angular part 
of the fermionic equation of motion will be solved numerically, and an exact numerical value 
for the angular eigenvalue will also be determined and employed in the calculation. 

The outline of our paper is as follows: we start, in Section 2, with the presentation 
of the theoretical framework for our analysis, and the derivation of the radial and angular 
part of the spin- 1/2 equation of motion on the brane. The expressions for the transmission 
coefficients and the various Hawking radiation emission rates, adopted for emission on 
the brane, are reviewed in Section 3. In Section 4, the numerical techniques used to 
integrate the angular and radial parts of the spin- 1/2 equation of motion and determine 
the angular eigenvalue are described. Section 5 presents our numerical results for the 



-3- 



transmission coefficient, the particle, energy and angular momentum fluxes, as well as the 
total emissivities, and their dependence on both the dimension of spacetime and angular 
momentum of the black hole is investigated. The angular distribution of the particle and 
energy emission rates are also studied by using the exact values of the spin-l/2-weighted 
spheroidal harmonics. We finish with a summary of our results and conclusions, in Section 
6. 



2. Fermion field in a brane-induced rotating black hole background 

We start our analysis by presenting the line-element that describes the 4-dimensional grav- 
itational background on a brane placed in the vicinity of a higher-dimensional rotating, 
uncharged black hole. The latter (4 + n)-dimensional spacetime is given by the well-known 
Myers-Perry solution pC|] . The induced-on-the-brane gravitational background follows by 
projecting the higher-dimensional one on the brane, by fixing the values of all additional 
azimuthal coordinates to ir/2. Then, the brane line-element assumes the form 1 |7| 



fi \ 2 2afi sin 2 6 S 



2 2 , a 2 // sin 2 en 2 2 



(2.1) 



Xd6 z - ( r z + a z + ^ ) sin 2 e?(V 



where 



A = r 2 + a 2 --^ T , S = r 2 + a 2 cos 2 9. (2.2) 



The mass and angular momentum (transverse to the rt^-plane) of the black hole are then 
given by 

(n + 2)A n+2 2 

Mbh = n M ' J = — r~n M BH a , (2.3) 

lo7rG n + 2 

with G being the (4 + ra)-dimensional Newton's constant, and A n+ 2 the area of an (n + 2)- 

dimensional unit sphere, given by 

27r (n+3)/2 

An+2 = r[(n + 3)/2] • {2A) 

The 4-dimensional background ( |2,1| ) is the one 'felt' by all Standard Model fields, includ- 
ing fermions, which are restricted to live on the brane. As the reader notes, the higher- 
dimensional black hole has been assumed to have only one angular momentum component, 
instead of the [(n + 3)/2] possible ones. The reason for this is that, under the assumption 
that the brane is infinitely-thin, the colliding particles have a non-vanishing impact pa- 
rameter only along the usual three spacelike coordinates; therefore, the resulting angular 



In order for the black hole to be considered a classical object, its mass Mbh will be assumed to be 
at least a few times larger than the fundamental scale of gravity M*. Then, under the assumption that 
the brane self-energy is of the order of M*, its effect on the gravitational background can be considered 
negligible. Analyses that address the complete bulk-brane-black-hole system have only led, up to now, to 
consistent brane trajectories solutions |n| - for an attempt to analytically construct a small black hole 
stuck on a tensional brane, see [B2l. 
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momentum will be around an axis also in our brane. An upper bound can be imposed 
on this sole angular momentum parameter of the black hole by demanding the creation of 
the black hole from the collision of the two particles. The maximum value of the impact 
parameter between the two particles that can lead to the creation of a black hole was found 
to be § 

'n + 2" 



1 + 



This, in conjunction with the relation J 
black hole, leads to 



(n+l) 



1 

(n + l) 



fi^ + L) . (2.5) 

bMsH /2 |22j for the angular momentum of the 
= (2.6) 



where a* = a/r^. For n > 1, the black hole horizon is given by the unique solution of the 
equation A(r) = 0, and may be conveniently written as r^ +1 = p/(l + a 2 ). 

Using the Newman-Penrose formalism, the Kinnersley tetrad and Boyer-Lindquist co- 
ordinates, Teukolsky [|33| decoupled the field equations for spin s = 0,1/2,1 and 2 field 
perturbations of a Petrov Type D (vacuum) background metric (that also includes the 
traditional 4-dimensional Kerr spacetime). Teukolsky expressed all the spin-field equations 
as one single 'master' equation, with the 'helicity', or spin-weight, h = (+\s\,—\s\) of 
the field as a parameter. As Teukolsky 's original equation does not apply to the case of 
propagation in a projected-on-the-brane higher-dimensional background, as the one given 
in Eq. (|2.l| ), here we have re-derived this equation. The new 'brane master' equation now 
describes spin s = 0, 1/2 and 1 field perturbations of both 'helicities', and has the form 



(r 2 + a 
A~ 



2\2 



2 ■ 2 

a sin f 



8 2 Q h 2ap d 2 n h 
dt 2 + Ar™- 1 dt 86 



+ 



1 



A~ h ^ | A 







2h 



dr 



r + p 
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sin ( 
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2h 



sin' 

aA'(r) 
2A 



8 2 n h 



+ 



l COS I 



sin 



dt 



- + h [/icot 2 9 - 1 + (2 - A"(r))6, 



h,\hl 



+ 



(2.7) 



X-T h , 



where p = r+ia cos 9. In the above, $7^ = 0^(4, r, 9, 6) represents the spin-field perturbation 
and Th the source term in the field equations - we refer the reader to 0, |33| for details 
on the derivation of the above equation, and for the precise definitions of ^(i, r, 9, 6) and 
T/j, in each spin case. Equation (|2.7| ) is formally the same as Teukolsky 's original four- 
dimensional equation, except for the different expression of the metric function A(r) which 
is now n-dependent, and the term with A"(r) which is identically zero in four dimensions 
(this was already noted in p3] for the spin-1 case) and, obviously, for h < 0. 

The 'master' equation ( |2.7| ) is clearly separable for any value of the 'helicity' h, so that 
its solution can be written as a sum over the Fourier modes: 



/+oo +°° +1 
duj ^2 ^2 h a Ah n A (t,r,9 
-OO , hi i 



l=\h\ m=—l 



(2.8) 



h n A (t,r,9,6) = hRA(r) h S\(9)e 



—iuit +imd 
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where h a A are the Fourier coefficients, and the set of 'quantum' numbers is denoted by 
A = {Imw}. 

The radial and angular ODEs resulting from the separation of variables of the 'master' 
equation are, respectively: 

~K 2 - ihKA'(r) 



A' 



d ( Ah+l d h R A\ 



dr 



A 



dr J 



+ 



A 



+ Aifwjr + h(A"(r) - 2)5 hj \ h \ - h X A 



where we have used the definition 



K = (r 2 + a 2 )uj — am, 



(2.9) 
(2.10) 



and 



d_ 

dx 



dx 



■i" )— I + a 2 ui 2 {x 



1) — 2haujx 
m + hx) 2 



1 



+ h^A + 2mauj + h 



h S^{x) = 0. 



(2.11) 



where x = cos 9. In the above, ^Aa is the constant of separation between the angular and 
radial equations. 

In the spin- 1/2 case, the wave function of the fermion is represented by two spinors 
P A and Q A ' [34], that, in the massless case, satisfy the equations V aa'P A = V AA'Q A = 0. 
By using the Kinnersley tetrad (i^, n^, m^, m*) and the Fourier decomposition (|2.8|) , these 
equations can be separated into first-order ODEs with respect to the 9 and r coordinates. 
For the 2-component spinor P A = (P ,^ 1 ), we make the decomposition 



P° = (r-ia cos 9)- 1 _ 1/2 R A (r) _ 1/2 S A (0) e^e^ 
P 1 = +1/2 R A (r) +1/2 S A (9)e im *e-^. 

Defining new radial functions 

P+ = A 1 / 2 +1/2 R A , P_ = _ 1/2J R A , 

the first-order radial equations take the form 

A 1 / 2 (d r - iK/A) P_ = XP+ 
A 1 / 2 (dr + iK/A) P + = AP_ 



(2.12) 
(2.13) 

(2.14) 



(2.15) 



where A = \f -\/ 2 \a- It may easily be shown that similar definitions for the components 
of the Q A> spinor lead to the same system of first-order radial ODEs. The coupled radial 
equations (2.15) may then be written in matrix form as follows 



dP 

dr 



iK A 

aT^ + avs- 



<Tl P. 



(2.16) 



where P = f p + ^ > anci °"l = ( i o ) anc ^ °"3 = ( o -i ) are ^ ne P aun s pi n matrices. By taking 
the complex conjugate of this equation it is clear that, if P is a solution, then so is criP*. 
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This transformation corresponds to time-reversal. On causal grounds we will only admit 
one set of solutions: those that are ingoing at the horizon. This alternative, but equivalent 
to Eq. (p.Sj), set of radial equations will be used for the numerical determination of the 
radial part of the spinor wave function, as we will see in detail in Section 4. 



3. Hawking radiation in the form of fermion fields 

We are interested in the Hawking radiation of an evaporating black hole on the brane, which 
is described by the ( "past" ) Unruh vacuum | U~ ) . The quantum field theory of fermions on 

I, but for 



the Kerr black hole has been extensively studied in the literature [35, 36, 37 



ease of reference, we now outline the key features. For more details of the general formalism 
relating to fermions on curved space-time, see [p9f , or [40 for a more recent review. 

In constructing the state \U~), it is most straightforward to use Dirac 4-spinors, in 
which case the wavefunction of the massless fermionic field (that, for historical reasons, we 
call 'neutrino') takes the form [35, |3(1, 37 j 



g—iuitgirrvp 




where 



. sin 2 9 p 2 } * 



Here L is the lepton number 
neutrinos. The two-spinor rj takes the form 35, 36 



Ri(r)Si(0) 
R 2 (r)S 2 (e) 



T= [A; 

||, so that L = +1 for neutrinos and L 



(3-1) 

(3.2) 
•1 for anti- 

(3.3) 



The angular functions Si and S 2 satisfy the equations |35|, 3f| |37|] 



d 

— + ( auj sin ( 
a9 



auj sin ( 



m 
sin 6 

m 
sin 6 



Si(0) = \S 2 (6), 
S 2 (9) = -XSi(9), 



and are related to the spin- weighted spheroidal harmonics hS\ ( 2.11 ) by j35] 



Si(0) 



sm( 



S A (co S 9), S 2 (9) 



8111 ( 



S*a(cos 9). 



(3.4) 
(3.5) 



(3.6) 



There is a subtlety in the definition of the radial functions Ri, R 2 , depending on whether 
we are considering neutrinos or anti-neutrinos. For neutrinos (L = +1), we have j 



3§|7] 



R^r) = P- 



R 2 (r) = P + 



but for anti-neutrinos (L 



•1), the radial functions swap over [37]: 
R l (r) = P+, R 2 (r) = P_. 



(3.7) 



(3.8) 
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This is particularly important when calculating the particle fluxes for neutrinos and anti- 
neutrinos. 

The usual "in" and "up" modes are defined by the behaviour of the radial functions 
near the horizon and at infinity [35, 36 1: 

1 



1 

/IP 



r 
r 



(Bl p e iujr * , 0) 



',A^e~ iQr *) r 



cc, 
rh, 

► OO, 

> r h , 



(3.9) 
(3.10) 



and correspond to modes originating either at past infinity or past horizon, respectively. 
In the above, r* is the standard "tortoise" co-ordinate, defined by 



and 



with 0, the horizon angular velocity 



dr 



UJ 



n 



r 2 + a 2 
AW ' 



u — mQ,, 



a* 



(3.11) 
(3.12) 

(3.13) 



{l + a 2 )r h 

In Eqs. (|3!9|) -( ^To|) , A™, B™, A^, B^ are constants of integration which satisfy the 
Wronskian relations [36] 



1 - \A\ 



\b: 



Al ! 



R up 



rjin 



(3.14) 



where • is either "in" or "up" . These relations are useful for simplifying the flux formulae 
we shall derive in this section. 

The ("past") Unruh vacuum \ U~) is defined in the usual way, by taking the "in" modes 
to have positive frequency with respect to co-ordinate time near infinity (which corresponds 
to to > 0) and using a set of modes having positive frequency with respect to Kruskal time 
in the vicinity of the past event horizon. Using the lemma in Appendix H of [41], a suitable 
set of modes having positive frequency with respect to Kruskal time near the past horizon 
is 



2 cosh 



uJ 



2T H 



cxp 



uJ 



4T H 



^A P + exp 



4T H 



\ ^down* 



(3.15) 



for any value of Co. The modes ^ own are given by the complex conjugate of the modes 
t/>a\ after changing the sign of the Kruskal co-ordinates. Similarly, a suitable set of modes 
having negative frequency with respect to Kruskal time near the past horizon is 

2 cosh | —7— j <, exp 



2T H 



-UJ 



4T H 



V>a P + exp 



4T H 



^down* 



(3.16) 
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for any value of u. The "down" modes are vanishing on the right-hand quadrant of the 
conformal diagram of the space-time, which is the region in which we are interested. There- 
fore we do not consider them further. Further details of this construction can be found in 
1 35, 41, 42]. We therefore find that a suitable expansion for the neutrino field in terms of 
our positive and negative frequency modes is (cf. J35"|): 




dw 



+ 



duo 



all L 



u)>0 



2 cosh 



2T H ) 



exp 



4T H J 



a A F + exp 



4T H J 



;upt 




(3.17) 



where —A = {£, — m, —uj} and Th is the Hawking temperature of the (4 + n)-dimensional, 
rotating black hole, given by 



(n + 1) + (n - 1) al 
4vr (1 + al) r h 



(3.18) 



The Fourier coefficients ^oa in (2.8) have been promoted to quantum operators d\ (where 
we have dropped the helicity subindex h for clarity). The state \U~) is then annihilated 
by the operators a\ and 6 A : 



al\U-) 



.up 



V = bf \u 



U' 



(3.19) 



We now have the machinery in place to compute the fluxes of particles, energy and angular 
momentum from the black hole. 

Firstly, the particle flux is given by [E! 



I de dip r 2 sin e(u-\r\u-}, 

at J a 

•J tJ CO 



(3.20) 



where J r is the radial component of the conserved current 



j» = -[*,>fy]. 

Here, the Dirac adjoint ^ is defined by [36, 39|: 



(3.21) 



(3.22) 



where 7 is the flat-space Dirac matrix. A suitable representation of the curved-space Dirac 
matrices 7^ for the space-time metric ( |2.ip can be found in |^6[ . We compute the fluxes for 
neutrinos v and antineutrinos v , using the field expansion ( p. 17 ) and the mode functions 



of the form ( |3.1| ). This gives the total particle flux per unit time to be [35]: 

Ta (-1/25*1 + +1/2-S 1 !) 



d 2 N v 



d(cos 6)dt 4tt 



j 00 1 

E E 



£=1/2 m=-t 



du 



oj=0 



exp (w/Th) + 1 



(3.23) 
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The transmission coefficient Ta is the transmitted proportion of the flux which, in the case 
of an "up" mode ( 3.10| ) is equal to 



T A = \Bl P \ . (3.24) 
By virtue of the Wronskian relations ( |3.14 ) this is the same as the absorption probability 



of an "in" mode. Similarly, the net particle flux is [35] 

= i % 2 ^ exp(^T H ) + l TA ~ ^ ■ (3 - 25) 

It is tempting, in view of Eqs. dJ|) and (ggD to view the contribution from neutrinos as 
being proportional to _ 1/2 5a an d that from anti-neutrinos as being proportional to +i/2<Sa- 
In fact, at infinity, this is the case, but only because of the form of the "in" and "up" modes 
in this limit and Eqs. (|3.7|) -( |3T8|) . It is not true more generally. 

The fluxes of energy and angular momentum per unit time are given in terms of 
components of the stress-energy tensor: 

^-=f d9 dtp r 2 sm6 (U-\T tr \U-), (3.26) 
f d9 dip r 2 sm6 {U-\T^\U~) , (3.27) 

J Srvn 



dJ 

~dt 



where the stress-energy tensor of the neutrino field is given by p 



T»u = - A ( [*,70»V„)¥] - [V (/1 *,7„)tf] )• (3.28) 



4 

The spinor covariant derivative is defined by [36] 



= [dp - r M ) % (3.29) 

where the T ^ matrices can be taken to be traceless and satisfy p(| 

7% ~ r„ 7 " + 7^ = 0. (3.30) 

These can be straightforwardly computed using the formula (cf. but with different 

conventions) 

r M = -^[7 ff ,7%]- (3-31) 

The formulae for the matrices are sufficiently large that we do not reproduce them here. 

After computing the expectation value of the stress-energy tensor in the state | U~ ) , 
the total fluxes of particles, energy and angular momentum emitted by the black hole per 
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unit time and frequency, in the form of fermion modes, are given by [35] 



exp (w/Th) 


+ 1 






exp (w/r H ) 


+ 1 


m 




exp (cj/Tr) 


+ 1 



dp Nis-^i/ 1 % % -\ 1 , , 

-jmT- = ^ 2^ 2^ PY n to /Ttt^ -i- 1 ' (3 - 32) 

<=i/2m=-« 
d 2 i£ 1 ^ v - W 

= 2^ 2^ 2^ P ™ /Ttt^ 4- 1 ' ^ 3 - 33 ^ 

<=l/2m=-< 
<=l/2ro=-< 

The emission rates ( |3.32| )-( |3.34| ) describe the different fluxes emitted by the black hole over 
the whole solid angle SI 2 . = 4ir, and follow by performing an integration over the angular 
coordinates 9 and (p. If we go one step backwards, we may derive the angular distribution 
of the emitted radiation by displaying the exact dependence of the differential rates on the 
azimuthal angle 9, and write 

d 3 N u+j? 1 °° 1 1 

~ai a\A+A = T~ X X /~ \ i i Ta (-l/^A + +1/2'S i a) > (3.35) 

a(cos 9)atauj gi/ 2 — t ex P \ ' R J 



d 3 £ 



1 x C<J 

/, m ;/ ; = T~ X! X! 7--/T \ ! i Ta (-i/2 g A + +i/2<Sa) ■ (3-36) 

d(cos9)dtdtu ^ n e z^ 2 _ e ex P (W^h) + 1 

It is easy to check that the above two expressions reduce to Eqs. (|3.32|) - (|3.33| ) when we 
use the angular function normalization 

J |±i/2Sa(cos60| 2 d(cos9) = J \ ±1 / 2 Sa(x) | 2 dx = 1 . (3.37) 



Note that in the particle ( |3.35 ) and power ( |3.36| ) fluxes, both 'helicities' h = —1/2 and 



h = +1/2 make a contribution. As we will see, this will lead to radiation spectra that are 
invariant under the parity transformation 9 — > ir — 9 (for more discussions on the asymmetry 
under parity of the derived spectra in purely four dimensions, see, e.g. [35, 37, 38 1). 



4. Numerical analysis 

The transmitted flux in each angular mode may be calculated by numerically solving the 
radial equations. Before we may begin, however, we must first solve the angular equations 
to determine the eigenvalue A = y / _ 1 / 2 AA that appears in the radial equations. In order 
to compute the angular distributions of the particle and energy fluxes, we will also need 
the angular eigenvectors: the spin-l/2-weighted spheroidal harmonics. In this section, we 
briefly outline the numerical methods used to solve both the angular and radial equations, 
and describe our method for determining the transmission coefficient Ta- 

The angular eigenvalues ^Aa were calculated using a spectral decomposition method, 
described in a paper by Hughes [44]. In appendix A of [44], it is shown that spheroidal 
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harmonics may be expanded as a sum over spherical harmonics of the same spin-weight, 



oo 



h S A (0)= ^ h b<ff h Y jm (9), (4.1) 

j — ^min 

where l m i n = max(\h\, \m\). A recurrence relation for the expansion coefficients is 
generated by substituting (4.1) into the angular equation ( |2.11| ) and exploiting the orthog- 



onality properties of the spin-weighted spherical harmonics. Writing bj = ^bf^ for clarity, 
the result is 

-h£jmwbj = [(aw) 2 cjl 2 j 2 ] bj-2 + [(aoj) 2 cf_ ltj7 2 - 2auhcf_ 1 j 1 ] bj-i 
+ [(auifc^ - lawhc^ - j(j + 1)] bj 

+ [(au) 2 cY +1>jj2 - 2auhc™ +ljj>1 ] b j+1 + [(au;) 2 cf +2J|2 ] b j+2 , (4.2) 

where we have defined 



< J j.i,2 = ^J^^{j t 2,m,0\l,m)(J,2,-h,0\l,-h) + ! L.:i) 



<%i,x= J^^(j,l,m,0\l,m)(j,l,-h,0\l,-h). (4.4) 

Here, m, 0\l, n) are Clebsch-Gordan coefficients. The separation constant h^Ai appear- 
ing on the left-hand-side of Eq. (|4.2| ), is related to the angular parameter h\\ by 



h £ A = h \ A + 2mauj - a 2 co 2 + h(h + 1). (4.5) 

In the spin-half case, £ and m are half-integers, and h = —1/2. In the non-rotating limit 
the solution reduces to Ei m = l{t + 1) and _i/ 2 Aa = (£+ |) 2 . 

The recurrence relations (f|^) may be rewritten in matrix form. The eigenvalue of the 
matrix is (minus) the separation constant h£\, and the eigenvector is a column of expansion 
coefficients bj. Since the matrix is band-diagonal, the eigenproblem may be solved in a 
numerically efficient way. For more details on this method, see Appendix A in p4|]. 

The results for the angular eigenvalues were independently checked using a partial 



fraction method [45], and a variant of the 'shooting' method [ }46[ , which is described in 
detail in |25|, ^7|. The angular eigenvalues coincided up to double precision. The values of 
the corresponding separation constants h^A following from the above analysis are shown 
in Fig. ||(a), for some indicative spin-1/2 modes. 

The spectral decomposition method may also be used to calculate the angular func- 



tions hSA(x), according to Eq. (O). However, this calculation would first require the 
determination of the spin-weighted spherical harmonics hYj m (9). Instead, we preferred to 
follow the series expansion approach introduced by Leaver ]45| . A solution to Eq. ( 2.11| ) 
may be expressed as 

l l 00 

h S A (x) = e auJX (l + x)2l' m - h l(l - x ) 2 |m+ft| Cn(l + x) n , (4.6) 

n=0 
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C0S(8) 



Figure 1: (a) The separation constant h^A, for the modes I — 1/2,3/2 and m = —£,...,£, and 
'heUcity' h = —1/2; (b) Spin-1/2 weighted spheroidal harmonics for a spinor field with h = —1/2, 
I = 1/2, m = -1/2, and for auj = (0, 1,2,3,4). 



where x = cos(#) and the expansion coefficients c n are related by a three-term recurrence 
relation, 



aoci + P c = 0, 

a n C n +l + PnC„. + 7n c n-l = °> 71 = 1,2,... 

The recurrence coefficients a n ,(3 n , 7„ are 



(4.7) 
(4.8) 



a n = -2(n + l)(n + 2A;_ + 1) 

f3 n = n (n - 1) + 2n(k- + k + + 1 - 2aw) 

- [2aw(2A;_ + h + 1) - + 
7 n = 2auo{n + k^. + k+ + h), 



+ k + + 1)] - [aV + fc£ A ] 



(4.9) 



with fe_ = ^l^n 



/i| and fc+ = !|m + This expansion, together with the eigenvalue h^A 
calculated via the spectral decomposition, provides an efficient method for the calculation 
of spin-half spheroidal harmonics, such as those shown in Fig. |l|(b). Note that _i/2>Sa i s 
clearly greater on the Southern hemisphere than on the Northern one, with a similar (but 
opposite) asymmetry arising also for +i/2S\- However, their sum is invariant under the 
change 9 — > tt — 8, which means that, unless both 'helicities' h = ±1/2 are taken into 
account in the calculation of the angular distribution of the radiation fluxes ( |3.35| ) and 



( 3.36 ), the resulting fermionic emission will be asymmetric with respect to parity. 

Once the eigenvalues of the angular equation are known, we may proceed to solve the 
radial equations. One possible approach is to start from the second-order radial Teukolsky 
equation ( |2.9[ ), and introduce a new radial "tortoise" coordinate (like the one in ([3.11 )) to 
transform to a potential barrier problem 

'' 2 V(r) 



R{r) = 0. 



(4.10) 



Equations of this form are familiar from quantum mechanical problems, and standard 
methods for finding transmission coefficients may be employed. However, a disadvantage 
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of this method is that the definition of the tortoise coordinate depends on the number of 
extra dimensions, and the relation r*(r) can be hard to invert (see the discussion in p5|1). 
Here, we prefer instead to work with the pair of coupled first-order equations ( p. 16 ). These 



may be solved by combining series expansions at the horizon and at infinity with numerical 
integration, as described below. 

Briefly, our method for determining the transmission coefficient Ta is as follows. We 
start by demanding that the solution be purely ingoing at the horizon. The ingoing solution 
may be written as a power series in (r/r^ — l) 1 / 2 (see next subsection). We evaluate this 
series solution at tenth order at a point close to the horizon (r/r/j — 1 < 10 ) and use this 
as the initial condition for a numerical integration scheme. We integrate outwards into the 
asymptotically-flat region, r/r^ > 100, using a fourth-order Runge-Kutta routine. Here, 
the solution is matched onto ingoing and outgoing power series solutions in r^/r which are 
also evaluated at tenth order. The asymptotic solution at infinity is a linear sum of ingoing 
and outgoing waves, 

A™ P in + A out P out , asr^+cxD, (4.11) 

where A m and ^4° ut are complex constants, and P in and p out = <riP in * are series expansions 
described below. Using our numerical solution we may determine A m and A ont up to an 
overall magnitude and phase. The transmission coefficient is then simply 



Ta = 1 



.4° 



A 1: 



(4.12) 



4.1 Expansion at the horizon 

We may write K and A as power series around r = r^, 

k=0 v h 7 k=l v h ' 

The first terms are 



(4.13) 



K = {r 2 h + a 2 )u - am, K x = 2ujr\, (4.14) 

A x = 2r\ + (r\ + a 2 )(n - 1), A 2 = r\ - \n{n - \)(r\ + a 2 ). (4.15) 

Defining rj = (r/r^ — l) 1 / 2 , the radial equations ( 2,16| ) may be written as 



y-MWP, (4.16) 
aij r] 

where 

oo 

Mfa) = M fc v h = (2wo) ^3 + V (7^) °i + ^ 2 2^o (g - sf ) °"3 + ■ • ■ (4.17) 

k=0 

and so = ifo^/i/Ai- This form makes it clear that r = is a regular singular point of the 
radial equations. 
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The ingoing solution at the horizon may then be written as 



P_ 
P + 



- 1 



-IS 



,/£7(l-4iso) 
1 + 



r? + ... 



h 



Ai(l-4is ) 



(4.18) 



Note that the P_ series expansion includes only odd powers of whereas the P+ expansion 
includes only even powers. The outgoing solution may be generated by time-reversal, that 
is, P out = o-iP 111 *. 

4.2 Expansion at infinity 

Similar techniques may be applied to find an expansion around the singular point at infinity. 
First, we define z = r^/r. The radial equations can then be expressed as 



dP 

dz 



p. 



(4.19) 



This makes it clear that the point at infinity is an irregular singular point (of rank 1) and 
we should look for a series expansion of the form 



V = e t ^ 1 z- t ^i k z\ 



(4.20) 



fc=0 



where to an d t\ are constants, and £ fc are two-component eigenvectors to be determined. 
Substituting this into the radial equations and equating the coefficients of like powers of z 
yields a system of equations 



(N + t )£o = 

(N + to)£i = -(Ni + ti)£ 

(N + 1 )£ 2 = "(Ni + (h - - N 2 ^ 



(4.21) 
(4.22) 
(4.23) 



k-2 



(N + t )£ k = -(N 1 + (h — (k— im k -i ~ E N *-^j 

3=0 



(4.24) 



The solution for to and t\ depends on the number of extra dimensions (n). In 4D, we find 

t = ±iu;r h , ti = ±iuj{r1 + a 2 )/r h , (4.25) 

with the positive choice of sign corresponding to an outgoing solution. In 4D such a series 
solution is given by 



pout 



e iujr r iu(rl+a 2 )/r h 



+ 



2tor 



A 2 - 2w 2 (r 2 + a 2 ) 2 /r 2 + 2auom 
-A 



+ 0[± (4.26) 
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In the asymptotic limit, the upper component dominates over the lower component in the 
outgoing solution. Again, the ingoing solution can be found by interchanging P_ and -P+ 
and taking the complex conjugate. 

In higher dimensions the index to is unchanged, but t\ = 0. The absence of a log- 
arithmic phase factor for higher-dimensional solutions indicates that the influence of the 
gravitational potential is now "short-range" for large r (up to r ~ L, where L is the size of 
the extra dimensions). The potential is still "long-range" near the horizon, however, and 
numerical troubles were avoided by using series expansions at the horizon (and infinity) up 
to tenth order. 



5. Numerical results 



In this section, we present exact numerical results for the transmission coefficient and 
the corresponding particle, energy and angular momentum emission rates for a higher- 
dimensional, rotating black hole emitting spinor fields on the brane. Their dependence 
on the dimension of spacetime n + 4 and the angular momentum parameter a* will be 
investigated, and the angular distribution of the emitted number of particles and energy 
will also be derived. In the final part of this section, the total emissivities of flux, power and 
angular momentum will be computed, and their dependence on the fundamental parameters 
n and a* will also be discussed. 

5.1 Transmission coefficient 

The first issue we need to address is the behaviour of the transmission coefficient Ta in 
terms of both the angular momentum parameter a* and the number of extra dimensions n. 
The value of Ta follows from the numerical integration of the radial equations ( p. 16 ) after 



the exact value of the eigenvalue ^Aa has been determined, as described in the previous 
section. 

The behaviour of the transmission coefficient also strongly depends on the particular 
mode studied, therefore in Figs. || and || we present the dependence of two indicative 
modes (£ = 1/2, m = 1/2) and (£ = 1/2, m = —1/2) on a* and n, respectively. Figure ||(a) 
reveals the monotonic behaviour of Ta for modes with m > that leads to a continuous 
enhancement in terms of the angular momentum parameter a*. On the other hand, as 
we may see in Fig. ||(b), modes with m < exhibit a non-monotonic behaviour that 
strongly depends on the energy regime that is probed: whereas the transmission coefficient 
is enhanced in the low-energy regime, it is suppressed in the high-energy one. The two 
panels of Fig. || depict the dependence of Ta on the number of the transverse-to-the-brane 
dimensions n. The exhibited behaviour is again found to be mode- and energy-dependent. 
Modes with m < have a transmission coefficient which is monotonically suppressed with 
n, while the behaviour of modes with m > is characterised by enhancement in the 
low-energy regime and suppression in the high-energy one. 

A similar behaviour in terms of a* and n was found for the corresponding modes of 



the emitted gauge bosons studied in [25]. The dependence of modes with m < is easily 



observed to be identical in the two cases. For modes with m > 0, though, the agreement 
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Figure 3: Transmission coefficients for spin-f/2 emission on the brane from a black hole with 
a, = 1 and variable n, for the modes (a) £ = 1/2, m = 1/2, and (b) £ = 1/2, m = —1/2. 



is not obvious since, in the case of gauge bosons, modes with m > are superradiant ones. 
Nevertheless, when the fact that a 'more negative' value of Ta in the superradiant case 
is equivalent to a 'more positive' value in the non-super radiant one, agreement is again 
found between the two cases. In [25], a detailed analysis was performed to show that 
the above behaviour of the transmission coefficient was directly linked to the behaviour of 
the effective potential seen by the propagating emitted mode - as the behaviour found for 
fermions is similar to those for gauge bosons, we will not repeat the same discussion here 
but, instead, refer the reader to Ref. [25| for more information on this point. 

The spectrum of the emitted fermionic modes will also be strongly determined by the 
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behaviour of the thermal factor (e^^ 11 + l)" 1 in terms of the same parameters n and a*. 
It may be easily seen that, for fixed a* and r^, the temperature of the black hole increases 
with the number of extra dimensions n, thus giving a boost to all emission spectra. The 
dependence on the angular momentum parameter a*, on the other hand, is much more 
subtle, and depends on the particular mode and energy regime. For modes with m < 0, a 
simple numerical analysis reveals that the exponent w/Tjj is an increasing function of a*, 
for all values of u, as long as n < 1, thus leading to a possible suppression of the emission 
rates for these modes as a* increases; however, for n > 1, this exponent is found to be a 
decreasing function of a*, for low enough values of the energy and high enough values of 
the angular momentum of the black hole. As the transmission coefficient is also enhanced 
in these regimes, according to Fig. |2](b), fermionic modes with negative m are more likely 
to be found in the low-energy part of the spectrum of rapidly rotating black holes. Finally, 
for modes with m > 0, the exponent is found to be a decreasing function of a* in the low- 
angular-momentum regime, and an increasing function otherwise, independently of the 
energy of the mode; as the transmission coefficient is also monotonically enhanced with a*, 
we expect these modes to mainly dominate the radiation spectra of slowly rotating black 
holes. 



5.2 Flux emission spectra 



In the next three subsections, we present the various emission spectra for spin-1/2 fields 
by using the definitions ( 3.32| )-( [3.36 ) and the exact numerical results for the transmission 
coefficient discussed above. We start with the particle flux, i.e. the number of fermions 
emitted by the black hole on the brane per unit time and unit frequency. Figures |] and 
|5| depict our results as a function of the dimensionless energy parameter ur^ and for fixed 
dimension of spacetime and angular momentum of the black hole, respectively. 




Figure 4: Flux emission spectra for spin-1/2 particles on the brane from a rotating black hole, for 
(a) n = 1, and (b) n = 6, and various values of a*. 



The dependence of the particle flux spectrum for fermions follows the same pattern as 
the one for either scalars |^, 24 1 or gauge bosons 2Sj. As either the angular momentum of 
the black hole or the dimension of spacetime increases, the number of particles emitted by 
the black hole is greatly enhanced. According to Fig. ||(a), for low values of n, the black 
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Figure 5: Flux emission spectra for spin-1/2 particles on the brane from a rotating black hole, 
with a* = 1, for various values of n. 



hole clearly tends to emit more low-energy particles than high-energy ones, although as a* 
increases the number of the latter becomes significant. For higher values of n, however, as 
Figs. |](b) and || reveal, the black hole consistently prefers to emit particles with energies 
wr^ > 1, with any increase of a* resulting in the shift of the emission curve towards even 
larger values of wr^. 

The angular distribution of the flux spectra of fermionic modes - characteristic in the 
case of a rotating decaying black hole - was also studied. It was found to borrow features 
from the corresponding angular distributions of both scalars and gauge bosons. The two 
indicative cases of n = 1 and n = 6, for three different values of the angular momentum 
parameter, i.e. a* = (0,0.6,1), are presented in Figs. ^ and [?]. The depicted spectra 
correctly reproduce the dependence on a* and n discussed above, and in addition reveal the 
existence of two forces that shape the final angular distribution pattern. The centrifugal 
force, acting on all particles independently of their spin, forces the emitted modes to 
concentrate on the equatorial plane, i.e. around the value 9 = tt/2 of the azimuthal angle. 
On the other hand, the non- vanishing spin of the particle leads to a spin-rotation interaction 
with the angular momentum of the black hole and tends to polarize radiation along the 
rotation axis (8 = 0,ir). The latter effect was, as expected, absent in the case of scalar 



fields [24] but quite prominent in the case of gauge bosons [25|. Here, as expected, it is 
also present but of a lesser magnitude due to the smaller value of the spin. According to 
Figs. |6] and [7|, and in agreement with the behaviour found for gauge bosons, the spin- 
rotation coupling is found to be particularly effective for particles with low energy, but its 
effect becomes subdominant compared to the one of the centrifugal force as either a* or n 
increases - the dying out of the polarization of the radiation along the rotation axis, as the 
angular momentum of the black hole increases, was also noted in the pure 4-dimensional 
case [35]. 



- 19 - 




Figure 6: Angular distribution of the flux spectra for emission of fermions on the brane from a 
rotating black hole, for n—1 and a* = (0, 0.6, 1). 




Figure 7: Angular distribution of the flux spectra for emission of fermions on the brane from a 
rotating black hole, for n = 6 and a* = (0, 0.6, 1). 

5.3 Power emission spectra 

We now turn to the power spectrum, i.e. the energy emitted by the black hole in the form 
of fermions per unit time and unit frequency. Our results are displayed in Figs. |8| and 
again for fixed dimension and angular momentum parameter, respectively. For fixed n, 
as the angular momentum of the black hole increases, the energy emission curve steadily 
moves towards larger values of the energy parameter due to the increasingly larger number 
of high-energy particles emitted, as was observed in the previous subsection. 

As a result, significantly larger amounts of energy are radiated away in the high-energy 
channels than in the low-energy ones compared to the case of a non-rotating black hole. At 
the same time, each curve becomes taller and wider, thus leading to a total energy emission 
rate that is substantially enhanced. According to Fig. ||, a similar, but significantly more 
important, enhancement is observed when the dimension of spacetime increases while the 
angular momentum is kept fixed: a mere comparison of the area outlined by the curves 
in the two extreme cases of n = 1 and n = 7 suffices to reveal the enhancement in the 
emission rate, as the number of the transverse-to-the-brane dimensions increases 2 . 

2 A comparison between our results and the ones produced in |^7| for the emission of fermions reveals a 
discrepancy of 3-5%. This might be due to either the failure of the analytic approximation to the angular 
eigenvalue used in [p7|, the potentially different number of partial modes included in the calculation, or 
simply the different methodology used to determine the transmission coefficient. 
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Figure 9: Power emission spectra for spin- 1/2 particles on the brane from a rotating black hole, 
with a» = 1, for various values of n. 



The angular distribution of the power spectra for fermion emission on the brane was 
investigated next, and the results are shown in Figs. [l(] and [Tl], for the indicative cases 
again of n = 1 and n = 6, respectively, and for a* = (0,0.6, 1). For a* = 0, the spectrum 
has no angular dependence as expected, however, as soon as the angular momentum of the 
black hole is turned on, the effects of the centrifugal force and the spin-rotation coupling 
become also apparent here. The latter factor features in the angular distribution pattern 
only in the low-energy regime and only in the case of low dimension (n = 1) - as discussed 
earlier, as either n or a* increases, the low-energy part of the spectrum becomes increasingly 
more suppressed compared to the high-energy one which results also in the disappearance 
of the polarization of the radiation along the rotation axis. As in the case of the flux 
spectra, the centrifugal force prevails in the high-energy emission channels, and the bulk 
of the energy is clearly emitted along the equatorial plane. 
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Figure 10: Angular distribution of the power spectra for emission of fermions on the brane from 
a rotating black hole, for n = 1 and a* = (0, 0.6, 1). 

Power Power Power 




Figure 11: Angular distribution of the power spectra for emission of fermions on the brane from 
a rotating black hole, for n — 6 and a* = (0, 0.6, 1). 

5.4 Angular momentum spectra 

We now proceed to study the rate of emission of the angular momentum of the black 
hole, and its dependence on a* and n. The amount of angular momentum radiated by the 
black hole per unit time and frequency, is shown in Figs. |l2|(a,b) for n = 1 and n = 6, 
respectively, and variable a*, and in Fig. |l^ for fixed a* and variable n. 




2 4 6 8 10 5 10 15 20 25 30 



to r h to r h 

Figure 12: Angular momentum spectra for spin-1/2 particles from a rotating black hole, for (a) 
n=l, and (b) n — 6, and various values of a*. 
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Figure 13: Angular momentum spectra for spin- 1/2 particles from a rotating black hole, for a* = 1 
and various values of n. 



As in the case of scalars and gauge bosons, the angular momentum spectrum strongly 
resembles the power spectrum. As a* increases, while n remains fixed, the emission curve 
shifts to the right thus allowing for the loss of the angular momentum of the black hole 
through the emission of higher energy particles. As either a* or n increases, the emission 
curve becomes again wider and taller, thus revealing the fact that rapidly rotating black 
holes, or black holes living in a spacetime with a large number of transverse dimensions 
lose their angular momentum much faster. 

5.5 Total emissivities 

We finally turn to the total emissivities, that describe the emission of either particles, energy 
or angular momentum per unit time, over the whole frequency band. These total rates 
follow by integrating the differential emission rates ( |3,32 )-(3.34) over the energy parameter 



u>r/j. Here, this integration has been performed from zero to a maximum value of the 
energy parameter, appropriately chosen to cover all the emission for the cases studied. For 
all three types of spectra, that maximum value was found to be ur^ = 30. As a result, the 
computed total emissivities are accurate, and can provide exact enhancement factors for 
all three types of spectra in terms of a* and n. 

In Figs. |l4|(a,b), we present the total emissivities of all three types of emission rates, 
particle flux, energy and angular momentum, as a function of the number of extra dimen- 
sions n and angular momentum parameter a*, respectively. The first histogram reveals the 
enhancement of all three rates, as n increases, while a* remains fixed (a* = 1). According 
to our numerical data, as n goes from 1 to 7, the total number of fermions emitted by the 
black hole per unit time increases by a factor of 35.9, the total energy by a factor of 98.7, 
and the total angular momentum by a factor of 61.5. The second histogram, Fig. |l4|(b), 
gives the total emissivities when the dimension of spacetime is kept fixed (n = 6) and a* 
varies. As a* increases from zero to unity, a significant enhancement takes place again, 
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Figure 14: Total emissivities for spin-1/2 particle emission on the brane from a rotating black 
hole as a function of (a) n, for a* = 1, and (b) a*, for n = 6. 



although of a smaller magnitude: the number of particles emitted per unit time is tripled, 
while the energy emissivity is enhanced by a factor of 5.7. 

Let us finally discuss the angular distribution of energy emitted by the black hole, 
in the form of fermions, as given by Eq. ( 3.36| ), but integrated again over the frequency 
parameter. The derived power emissivity is then a function only of the azimuthal angle 
6, and for fixed n and a*, is depicted in Figs. Il5|(a,b), respectively. Similarly to the 




cos(e) 



Figure 15: Power emissivity (up to ujrh = 6) for spin-1/2 emission on the brane from a rotating 
black hole as a function of cosf5, (a) for n = 6 and variable a*, and (b) for a* = 1 and variable n. 



non-integrated one studied in a previous section, the power emissivity is symmetric with 
respect to the equatorial plane (cos 9 = 0). The effect of the centrifugal force on the angular 
distribution of the radiation is also obvious here with the bulk of emitted energy being 
concentrated on the equatorial plane, for all values of the angular momentum parameter 
and dimension of spacetime. However, the expected polarisation along the rotation axis 
caused by the spin-rotation coupling, that gave rise to a twin-peak pattern in a similar 



plot for gauge bosons [25|, is not visible here. This is due to the fact that the low-energy 
regime, where the spin-rotation coupling is effective, is significantly suppressed when the 
total emissivities are computed; in addition, the magnitude of this effect for fermions was 
much smaller in the first place, compared to the ones for gauge bosons, due to the smaller 
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value of their spin. The shape of the fermion power emissivity curves appearing in Figs. 
15 are found to agree with the ones derived in the case of a 4-dimensional, rotating black 
hole ||. 



6. Conclusions 

In this work, we have completed the study of the Hawking radiation on the brane from a 
higher-dimensional rotating black hole by investigating the emission of fermionic modes. 
A comprehensive analysis was performed that led to the particle, power and angular mo- 
mentum emission rates, and shed light on their dependence on fundamental parameters of 
the theory such as the spacetime dimension and angular momentum of the black hole. In 
addition, the angular distribution of the emitted modes, in terms of number of particles and 
energy, was thoroughly studied. Our results are valid for arbitrary values of the energy of 
the emitted particles 3 and angular momentum of the black hole, and complement previous 



results on the emission of brane-localised scalars [£3|, 24] and gauge bosons pqj . 

As in the case of our previous analyses for spin-0 and spin-1 particles, the equation of 
motion for fermions propagating in the projected-on-the-brane gravitational background 
needed to be determined first, a task that was performed by employing the Newman- 
Penrose formalism. The derivation of exact results for the various emission rates demanded 
the numerical integration of both the angular and radial part of the equation of motion 
in exactly that order. The former calculation offered exact results for both the angular 
eigenf unctions, the spin-weighted spheroidal harmonics, and the angular eigenvalues, that 
also appear in the radial equation. The latter calculation led to exact results for the 
radial part of the wavefunction of the fermionic field, and eventually for the transmission 
coefficient. 

Our numerical results for the transmission probability as well as for the particle, energy 

and angular momentum fluxes were presented in Section |5|. The differential emission rates 

per unit time and frequency, integrated over the azimuthal angle 6, were addressed first, 

and spectra for different values of the dimension of spacetime and the angular momentum 

of the black hole were constructed as a function of the energy of the emitted particle. As 

in the case of scalars and gauge bosons, it was revealed that any increase in either n or 

a* significantly enhances all emission rates of the black hole. The enhancement is present 

in all frequency regimes, and significantly boosts the emission of higher energy particles 

compared to the case of a non-rotating black hole. Total emissivities, i.e. emission rates 

integrated over all frequency regimes, were also presented at the end of Section || these 

total emission rates confirmed the enhancement of all spectra, in terms of n and a*, and 

led to the derivation of the corresponding enhancement factors: as n increases, the various 

emissivities are enhanced by a factor that varies from 35 to almost 100, while, as a* grows, 

this factor is of the order of 3 or 6. 

3 Using the results in E3], we have investigated at what value of cor^ is the contribution to the spectrum 
from ultra-short distances (where a full theory of quantum gravity would be required) of the same order 
as the total spectrum itself. For the main cases plotted in this paper, such value of uifh is at the tail-end 
of the curve for each dominant mode, so that the main contribution to the spectrum is indeed from the 
semiclassical regime. 
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Throughout the study of the emission of various spin-s particles on the brane by a 
higher-dimensional rotating black hole, a significant amount of energy was invested in 
the calculation of the exact spin- weighted spheroidal harmonics, that would allow us to 
determine the exact angular distribution of the emitted particles and energy. Contrary 
to the case of a spherically-symmetric black hole, the various emission rates for a rotating 
black hole bear a signature of a characteristic direction in space which is the axis of rotation. 
It was found that the majority of particles emitted by the black hole are restricted to do 
so along the equatorial plane due to the strong centrifugal force - this conclusion was 
reinforced after the total energy emissivities as a function of the azimuthal angle were 
also computed. Nevertheless, the spin of the particle can significantly affect the angular 
distribution pattern at the low-energy regime due to the spin-rotation coupling with the 
angular momentum of the black hole. For comparison, in Figs. |l(](a,b,c) we present the 
angular distribution of the energy emission spectrum for a 6-dimensional (n = 2) rotating 
black hole emitting radiation on the brane, in the form of scalars [24], fermions and gauge 
bosons p3], respectively. From these, it becomes obvious that the angular variation in the 
pattern of the emitted radiation can be a distinctive signature of emission from a rotating 
black hole, but also of the spin and energy of the particles emitted. 




Figure 16: Angular distribution of the power spectra for (a) scalars, (b) fermions, and (c) gauge 
bosons on the brane from a 6-dimensional black hole with a* = 1. 

The array of graphs displayed in Fig. [l6| can shed light also on the question of the 
relative emissivities of different types of particles on the brane from a rotating black hole. 



In the absence of angular momentum, it was shown [15| that the preference of the black 



hole for the emission of scalar fields in four dimensions is replaced by the one for gauge 



bosons when a number of additional spacelike dimensions is introduced. In [25], we were 
able to show that a higher-dimensional rotating black hole also prefers to emit gauge 
bosons over scalars on the brane. With the results for fermions being now available, the 
complete picture has emerged. By merely comparing the vertical axes of the graphs in Fig. 
|j~6| , one may see that a higher-dimensional rotating black hole emits approximately the 
same amount of energy on the brane in the form of scalars and fermions - a more careful 
comparison of the emission rates integrated over the azimuthal angle, presented in Fig. |8|, 
with the corresponding ones in |24| reveals that the fermions slightly dominate over the 
scalars. Nevertheless, as Fig. |l~6|(c) reveals, a significantly larger amount of energy (by, at 



least, an order of magnitude) is spent by the black hole in the emission of gauge bosons 
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the above result holds for all values of n. 

The results presented in this work, addressing modifications of the properties of black 
holes submerged into a higher-dimensional brane- world spacetime, have an obvious theoret- 
ical interest. The potential detection of the Hawking radiation from a higher-dimensional 
black hole formed during a high-energy particle collision in a ground-based accelerator, in 
the near or far future, increases their significance as the corresponding spectra offer vital 
information about fundamental parameters of the system under study, such as the dimen- 
sion of spacetime and angular velocity of the black hole, as well as for the creation of the 
black hole itself. This information is obtained by looking only at the emission channel along 
our brane where ordinary Standard Model particles, as well as the potential observers, are 
restricted to live. With the investigation of the spin-down and Schwarzschild phases of the 
life of the black hole now completed, what remains to be done is a realistic study of the 
dynamical evolution of the produced black hole where more advanced aspects, such as the 
variation of the temperature of the black hole as the decay progresses or the properties and 
interactions of the Standard Model particles, are properly taken into account. This goal 
will be achieved by entering our results for the emission of individual scalar, fermionic and 
gauge bosonic degrees of freedom from a constant-mass black hole into the CHARYBDIS 
black hole event generator [49]. We hope to report results from this analysis soon. 
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